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Abstract 

Discretizations of the Euler top sharing the integrals of motion with the continu- 
ous time system are studied. Those of them which are also Poisson with respect to 
the invariant Poisson bracket of the Euler top are characterized. For all these Pois- 
son discretizations a solution in terms of elliptic functions is found, allowing a direct 
comparison with the continuous time case. We demonstrate that the Veselov-Moser 
discretization also belongs to our family, and apply our methods to this particular 
example. 
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1 Introduction 



The subject of integrable discretizations of integrable dynamical systems is slightly more than 
twenty years old. During the first decade of its existence several different approaches were 
proposed for discretizing soliton equations ||AL|| , [||], [pJM| , ||QNCV|| , NCW| , but neither 
of them dealt with integrable systems of the classical mechanics. The first examples in this 
important subarea seem to appear about ten years ago jV], JSJ. The Veselov's paper contained 
discrete time versions of such classical integrable systems as the Euler case of the rigid body 
motion and the Neumann system. The algebraic construction behind these examples was 

PV |. 



later elaborated in more detail in 



However, the nature of these exceptionally beautiful 
examples remains somewhat mysterious. They still resist to be included in any of the existing 
general frameworks for integrable discretizations (cf. | |DL'1| |). The origin of the Veselov- 



Moser's discretizations, their place in the corresponding continuous time hierarchies, their 
relations to other existing examples, and several further points remain to be clarified. The 
problem of comparing the explicit solutions found in [|V|] , |[M V|| with the classical continuous 
time counterparts was also left open (despite the fact that on the level of equations of motion 
the continuous limit is easy to perform). 

The present work may be considered as a sort of an extensive comment on 



MV 



in 

the part concerned with the Euler top. We do not close all the open problems mentioned 
above, but we do hope to bring some light into some of them. In particular, we find formulas 
for the solution of the discrete time Euler top in the form which makes the comparison with 
the continuous time case immediate. Moreover, the Veselov-Moser system turns out to be 
by no means the only reasonable discretization of the Euler top. We introduce a whole 
family of discretizations sharing integrals of motion with the continuous time system, and 
characterize those of them which share also the underlying invariant Poisson structure. In 
this context the Veselov-Moser system becomes just one particular case, and not the most 
simple one. 

flMVfl 



It has to be mentioned that, unlike 



our procedure for obtaining explicit solu- 



tions is rather pedestrian; it does not use such advanced tools as spectral theory of difference 
operators @ or Baker- Akhiezer functions [MV|| . Actually, our construction is based on the 
addition formulas for elliptic functions, and could be invented already by Jacobi. However, 
as a heuristic tool for finding it we used the Lax representation of the Euler top in su(2) 
with a spectral parameter on an elliptic curve, which was certainly unknown in the times of 
Euler and Jacobi. 
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2 The Euler top 



The famous Euler's equations describing the motion of a rigid body with the fixed center of 
mass (Euler top) read |G): 

Mi = Q-^)m 2 M 3 , 

M2 = (^-£)M 3 M 1; (1) 
M 3 = Q-^MiMs. 

Here M = (Mi, M2, M 3 ) T is the kinetic momentum vector in the coordinate system attached 
firmly to the body; the axes of this system coincide with the principal axes of inertia, and 
the numbers A, B, C > are the corresponding moments of inertia. 
The vector form of the equations (|XJ) is: 

M = M x n(M) , (2) 

where the vector of the angular velocity is introduced: 

n , . /Mi M 2 M 3 \ T , . 

The equations (|l|) are Hamiltonian with respect to the following Poisson bracket on 
M 3 (M): 

{Mi,M 2 } = M 3 , {M 2 ,M 3 } = M 1; {M 3 ,M 1 } = M 2 . (4) 
with the Hamilton function H(M) = E(M)/2, where 

Ml M| M 3 2 
E(M) = -± + + -± . (5) 

The bracket (|j) is degenerate and has one Casimir function 

M 2 (M) = Ml + Ml + M 3 2 . (6) 

Generic symplectic leaves of the bracket (§) are two-dimensional spheres M 2 (M) = const, 
hence each system Hamiltonian with respect to this bracket is integrable in the Liouville- 
Arnold sense; in particular, the Euler top is integrable. 

An explicit solution to the equations ([I]) can be given in terms of the Jacobi elliptic 
functions 0. Suppose for the sake of definiteness that 

A > B > C > , (7) 
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then 



1 E 1 
— < < — . 

A~ M 2 ~ C 



The formulas for the solution look different depending on whether this quantity is greater 
or smaller than 1/B: 



Mi(t) = a 



cn u(t — t ) 
dn vit — t ) 



M 2 (t) = bmu(t-t ), M 3 (t) = c- 



dn vit — to) 
cn v[t — to) 



< ~Es < ^! while the lower ones refer to the case A < -J? < -g. The module k of the 



Here and below in similar situations the upper expressions in curly brackets refer to the case 
^ < -Es < 5, while the lower onee 
elliptic functions in (S) is given by: 

' (A-5)(M 2 -C£) ' 
(5- C){AE-M 2 ) 

(B — C)(AE — M 2 ) 
{ (A-B)(M 2 -CE) J 

so that we have always < k 2 < 1. The coefficients a, b, c are defined by: 



k 



(9) 



2 M 2 — CE u2 

a = A ^c~^ h 



B 



B 



M 2 -CE \ 
B-C 

AE-M 2 
A-B 



, 2 = C 



AE-M 2 
A-C 



(10) 



Finally, the frequency v is given by: 



v 



' (B-C)(AE-M 2 ) 
ABC 

(A- B){M 2 - CE) 
ABC 



(11) 



In the case -M^ = -g the elliptic functions degenerate to the hyperbolic ones, and we have: 



M,(t) 



cosh v(t — to) 



M 2 {t) = b tanh v{t - t ) , M 3 (t) 



cosh u(t — 1 ) 



where 



^ 2 _A(i?-C) M2 b2 = M 2 



B(A-C) 



B(A - C) 
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and 



2 (A - B)(B - C) - 
AB 2 C 



In all three cases the numbers a,b,c,v satisfy the condition abcv > 0. In what follows, we 
shall denote by i/(M 2 ,E) the positive square root of the function given in (|TTD, and hence 
we must have abc > 0. 



3 Discretizations 

Our aim here is to find integrable discretizations of the Euler top. 

Definition 1 An integrable discretization of the Euler top (fj) is a one-parameter family of 
diffeomorphisms F : M 3 x [0, e) i— ► R 3 , 

M = F(M,h), (12) 

such that each one of them F(-,h) : M 3 i— > M 3 is Poisson with respect to the bracket (j^j, has 
two integrals of motion M 2 and E, i.e. 

M 2 (F(M, h)) = M 2 (M) , E(F(M, h)) = E(M) , (13) 

and the following asymtotics hold: 

F(M,h) = M + hMxfl(M) + o(h), h^O. (14) 

Actually, our maps always will be defined by implicit equations of motion: 

M - M = hf(M, M, h) . (15) 

Proposition 1 Let f : M 3 x R 3 x [0, e) t— > M 3 be a C 1 -function in each argument such that 

/(M,M,0) = M x n(M) . 
Then for h small enough the equation ( |i~5|) defines a local diffeomorphism (|1^) satisfying 



Proof. This follows immediately from the implicit function theorem for all h satisfying 
det(J — hdf '/<9M) ^ 0, i.e. for all h small enough. ■ 

When speaking about discretizations, it is natural to think about M in as about 



sequences M m : Z i— > M approximating solutions M(t) : K i— > M of the Euler top (g) in the 
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sense that M m m M(m/i). In this context the formula (|i~2f) takes the form of the difference 
equation 

M m+1 = F(M m , h) , 

and, similarly, the formula ([15]) has to be thought of as 

M m+1 - M m = hf{M m , M m+1 , h) . 

The approximation property M m = M(mJi) + 0(/i) holds on finite time intervals and is 
assured by ([TJ|) (this is a standard fact from the numerical analysis). 



4 Elliptic coordinates and Poisson discretizations 



Now we would like to find a manageable criterium for a map fll2|) to be Poisson with respect 
to the bracket (^). The corresponding statement becomes rather transparent in a new 
coordinate system in R 3 (M). The corresponding change of variables 

$ : (M 2 ,E,p) ^ (M 1 ,M 2 ,M 3 ) 

is suggested by the formulas ([8]) for a solution of the Euler top. In the above formula 
(Mi, M 2 , M 3 ) G M 3 , and (M 2 , £7, cp) G M + x K + x T, where T = R/(4ifZ), and if = K(fc 2 ) is 
the full elliptic integral of the first kind corresponding to the value of k 2 given in (0). The 
formulas for $ are different in two regions of 1R 3 separated by the two planes 

E{M) - 1 « a-\)*=a-i)#. u«> 



M 2 (M) B \B A J \C B, 

Each one of the sets j A f i^j-, > ^} and | A f 1^ < ^} consists of two sectors, and each one 
of these sectors is bounded by two half-planes. 

The elliptic coordinates in R 3 (M) are introduced by: 

M 1= alTZ), M 2 = f,sn P , M3 = c|l^l. (17) 



dmp 

with a, 6, c given by (0), and the modulus fc 2 of the elliptic functions given by ([|). On 
the subset j-J^ > -^j the sign of c coincides with the sign of M 3 (which is constant in each 
one of two sectors), and the signs of a, b satisfy the condition sign(afe) = sign(c). Similarly, 
on the subset {^^y < ^} the sign of a coincides with the sign of M 1; while the signs of 
b, c satisfy sign(fec) = sign(a). Finally, on the four half-planes described by ([H]) the elliptic 
coordinates are defined by continuity according to 

M 1 = a , M 2 = b tanh (p , M 3 = — °- — , 
cosh (p cosh if 

the signs of a, c being the same as the signs of Mi, M 3 , respectively, and sign(6) = sign(ac). 
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Proposition 2 The formulas (|17j) define a valid change of variables (a local diffeomorphism) 
near each point o/R 3 . 

Proof. It can be verified by a direct calculation that in the region | jfl^ 7^ ^} we have: 

d(M u M 2 ,M 3 ) 1 



d(M 2 ,E,p) " Au(M 2 ,E) 



dn> ^ 



and that the partial derivatives of which this Jacobian is composed allow a continuation to 
the boundary {^f^gy = ^}. ■ 

Now it is obvious that an arbitrary map (0) having two integrals of motion M 2 (M) and 
-E(M), can be cast, in the elliptic coordinates (M 2 , E, (p), in the form 

M- M J . E = E, p = <p(M 2 ,E,p) . (18) 



Proposition 3 The map having two integrals of motion M 2 (M) and E(Wl) is Poisson 
with respect to the bracket (|4j) iff in the elliptic coordinates (M 2 ,E,p) it takes the form 

M 2 = M 2 , E = E , (p = p + g(M 2 ,E) (19) 

with the function g not depending on p. 

Proof. To prove this statement, we have to calculate the Poisson bracket (||) in the coordi- 
nates (M 2 ,E,ip). The corresponding formulas read: 

{M 2 , E} = {M 2 , <p} = , {E, if} = -2z/(M 2 , E) . (20) 

Indeed, the function M 2 is a Casimir function, hence it Poisson commutes with both E 
and if. To calculate the bracket {E,Lp}, we substitute (away from the boundary -M^ — ^) 
the expressions ( PD into an arbitrary one of the formulas (^), and after straightforward 
calculations arrive at the expression given above. 

Actually, the concrete expression for {E, if} is not essential for the proof of our proposi- 
tion. The only important thing is that this Poisson bracket does not depend on (p. Indeed, 
we have: 

= {EM- 

Since the Poisson bracket {E, ip} depends only on M 2 , E which are integrals of motion, we 
see that the necessary and sufficient condition for our map to be Poisson reads dip /dp = 1, 
which is equivalent to the last equation in (|T^). ■ 
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5 Special discretizations 

We derive now a family of discretizations of the Euler top. Our derivation is based on a Lax 
representation with a spectral parameter for the Euler top. We prefere to work with a su(2) 
Lax representation, since our experience in discretizing various geometric structures (see, for 
example, |BH| ) conviced us that this procedure may be performed most straightforwardly 
when applied to su(2) Lax formulations. 

To find a suitable Lax representation for the Euler top, we use a stationary version of 
the Lax representation of the chiral field model due to Cherednik [|Ch|| . Set 

I 3 

M ( u ) = ^E M k w k {u)a k , (21) 

II k=i 

where a k (k = 1,2, 3) are the Pauli matrices, and w k (u) are the following elliptic functions: 

/ \ 1 i \ dn(w, k) cn(w, k) 

wi{u) = p— r, w 2 {u)=p— -, w 3 {u) = p— -. (22) 

sn(w, K) sn(u, k) sn{u, k) 

Here the parameter p and the module k of the elliptic functions are defined by 

2 _ C(A-B) 
~ AC ' " B(A-C) ■ {2S) 

Further, for a vector V = (Vi, V 2 , V 3 ) T G R 3 set: 

I 3 

V(u) = — 22 V k w k (u - u )a k , (24) 

II k=i 

where the point u is chosen so that 

Wl (u ) = A' 1 ' 2 , w 2 (u ) = B- 1 / 2 , w 3 (u ) = C' 1 ' 2 . (25) 
Consider now the Lax equation 

M(u) = [M(u),V(u)] . (26) 

With the help of identities 

Wj{u)w k {u - m ) = Wj(u )wi(u - u ) - w k (u )wi(u) , (27) 

where (j, k, I) is a permutation of (1, 2, 3), one sees that the matrix equation (|26|) is equivalent 
to the set of the following two equations: 

M = M x n 1/2 (V) , (28) 
= V x ft 1/2 (M) . (29) 
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Here the second equation is equivalent to V = 7fi 1 / 2 (M) with some 7 G M, which, being 
substituted in the first equation, results in 

M = 7M x O(M) . (30) 

Obviously, for a nonvanishing function 7 = j(t) the latter equation is nothing but a time 
reparametrization of the Euler top (0). In this sense (|26| ) is a Lax representation of the Euler 
top. 

Discretizing the above construction in time, we introduce the matrix 

h 3 

V(u) = I + — 22V k w k (u-Uo)a k , (31) 
z% k=i 

and consider instead of ( ^6|) the discrete time Lax equation 

M(u) = V~ l (u)M(u)V(u) . (32) 

(this equation may be interpreted as a stationary version of the lattice chiral field model by 
JPfl). Representing the latter equation as V(u)M(u) = M(u)V(u) and using the identities 



(pTD, we find that our matrix equation is equivalent to the set of the following three equations: 

M-M = ^(M + M)x0 1/2 (V), (33) 
= V x n 1/2 (M + M) , (34) 

and 

3 

J2(M k - M k )V k w k {u)w k (u - u ) =0. (35) 

k=l 

Now (|34]) is equivalent to 

V = ^ 7 r2 1/2 (M + M) (36) 
with some 7 e R, which, being substituted in (|33|) , results in 

M-M = ^ 7 (M + M) x fi(M + M) . (37) 

It remains to notice that, plugging (|36|) , (^7|) into (|35|) , we bring the latter equation to the 
form 



3 

w k (u)w k (u - u )w k (u ) (wf(u ) - w|(«o)) = 

fc=l 
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which is automatically satisfied due to the identity 

Wi(u)wi(u - U )Wj(u ) - Wj(u)Wj{u - U )wi(u ) = W k (u )(wi(u ) - Wj(Uo)) . (38) 



So, the matrix equation (|32|) is equivalent to ( 37f) with some 7 G K, the vector V being given 

by®. 



However, the equation (]37|) does not completely define the discretization due to arbi- 
trariness of 7. In what follows we shall consider this equation with 7 being a certain fixed 
function on M, M. In other words, the subject of our further investigations will consist of 
discretizations governed by implicit equations of motion of the following special form: 

M — M = ~hj(M, M, h) (M + M) x fi(M + M) . (39) 

We call them special discretizations. It is obvious that if 

7 : R 3 x R 3 x [0, e) h-> R + 

is a C 1 -function of each argument satisfying 

7 (M,M,0) = 1, (40) 

then the equation (|39| ) fulfills the conditions of Proposition [I], and therefore defines a map 
(0)- 

In components, special discretizations may be presented as: 

M 1 -M 1 = \h-y(M, M, h) ( 4 - ^) (M 2 + M 2 ) (M 3 + M 3 ) , (41) 



4 ' v ' ' 7 Vfi C 
Mi Mi = i/i7(M,M,/i) [L-L]iM, + M,)iy[ v + M v ) . t 12) 

M 3 - M 3 = i/i7(M, M, /i) f 1 - ^ ) (Mi + Mi)(M 2 + M 2 ) . (43) 



4 ' v ' ' ' \A B 

6 Poisson property of special discretizations 

We investigate now the integrability properties of special discretizations (which are naturally 
expected due to the existence of a Lax representation with a spectral parameter). 

Proposition 4 Maps defined by equations of motion ( j3lf) possess M 2 (M) and E(M) as 
integrals of motion. 
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Proof. It follows from ( p9|) that 

(M — M, M + M) = and (M — M, f2(M + M)) = , (44) 

or, equivalently, 

(M,M) = (M,M) and (M, O(M)) = (M, O(M)) . (45) 

This proves our statement. ■ 

Actually, this statement can be almost inverted. Namely, consider an arbitrary discretiza- 
tion having M 2 (M) and E(N1) as integrals of motion. Then the pairs (M, M) satisfy fl4l)|). 
which is equivalent to ([44]). If, in addition, (M + M) x f2(M + M) ^ 0, then there exists a 
real number 7 = 7(M, M) such that (^) holds. 

In what follows we shall need also the fo lowing technical and obvious statement. 



Lemma 1 The set of fixed points of the maps ( p~2|) defined by (|39|) coincides with the set of 
the points M for which (M + M) x f2(M + M) = ; and this coincides with the union of the 
coordinate axes 

{M 1 = M 2 = 0} U {M 2 = M 3 = 0} U {M 3 = M x = 0} . 

For all other points M ; at least two of the three expressions Mi + M\, M2 + M2, M3 + M3 
do not vanish. 

Now we find conditions for a map defined by (|39D to be Poisson. To this end notice that 
the pairs (M, M) = (M, F(M, hfj belong to the subset of R 3 (M) x R 3 (M) singled out by 

the conditions M 2 (M) = M 2 (M) and E(M) = E(M). The elements of this subset may be 
parametrized by the quadruples (M 2 ,E,{p,(p) according to 

M=(.{^},^,c{^}), («) 
Equivalently, we can use the following coordinates: (M 2 , E, <p, A<p), where 

9 = ^-, *r = ^. (48) 

Let us denote on the above mentioned subset: 

7 (M,M,/i) = r(M 2 ,E,<^, A<p,h) . (49) 
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Theorem 1 The equations of motion of the special discretization (|39|) in the elliptic coordi- 
nates (M 2 , E, (p) have the form (|T8"D ; where the function (p(M 2 , E, ip, h) is implicitly defined 
by the equation 

T(M 2 ,E,y, A<p,h) = _J1 sn(Ay?) 

1 - Psn 2 (<^)sn 2 (Ay?) hu(M 2 ,E) cn(Ay?)dn(Ay?) [ ' 

The map defined by ( p9|) is Poisson, iff the equation fl5"0|) may be solved for Aip as 

2 

with the function S not depending on (p. 
Proof. Denoting 

D = D(M 2 , E, <p, Aip) = 1 - k 2 sn 2 (ip)sn 2 (Aip) , (52) 
we find from fliSD , ( f47j ) with the help of addition formulae for elliptic functions: 



Aip = 6(M 2 , E, h) = hu ( M ^ E ) (l + 0(h)) (51) 



2a j sn(^)dn(^)sn(A<^)dn(A^) 1 — _ 2a J cn(y?)cn(A^) 

D I fc 2 sn(^)cn((^)sn(A( / 9)cn(A( / 9) j ' 1 + 1 D 1 dn(^)dn(A^) 



(53) 

M 2 -M 2 = — cn(<^)dn(^)sn(Ay?) , M 2 + M 2 = — sn(^)cn(A^)dn(A</?) , (54) 

ft _ M k 2 sn(<p)cn(<p)sn(A<p)cn(A<p) \ ~ _ 2c f dn(^)dn(A<^) 1 

3 3 ~ D 1 sn(y?)dn(^)sn(A</?)dn(Ay?) J ' 3 + 3 £ \ cn(^)cn(A^) J ' 

(55) 

Plugging this into an arbitrary one of the equations of motion (f4T|)- (|43|), we arrive after 
some cancellations at the equation (|50D . Notice that Lemma [l] assures that, away from the 
fixed points, these cancellations are legitime in at least one of the equations of motion. 
So, we have arrived at the equation (^0]) of the form 

$>(M 2 ,E,ip, A(p,h) = , (56) 



which serves for determining the function (p(M 2 ,E,ip,h) for our map. The equation ((& 
may be rewritten as 

*(M 2 , E, ?±Z, h) d = f *(M 2 , E, <p,ip,h) = 0. (57) 

Proposition |3] gives a necessary and sufficient condition d(p / dip = 1 for our map to be Poisson, 
and this is equivalent to the following condition: 

-— + — — = on the solutions of (57) . (58) 
dip dip 
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Since, obviously, 

<9vj> _ 1 (dV dm \ dm _ 1 (dm dm \ 

dfi~2 \d$ + d(A<p) J ' Ify ~ 2 \ d$ ~ d(A<p) J ' 

we find that the Poisson property is equivalent to the following condition: 

dm 

— - =0 on the solutions of (56) . (59) 

In turn, ( |59| ) assures that the solutions of ( [56"] ) for Aip do not depend on (p. ■ 

Corollary. In the conditions of the above theorem, solutions M m = (Mi >m , M 2jm , M 3 m ) T 
to the difference equation 

M m+1 -M m = -h^(M m , M m+1 , h) (M m + M m+1 ) x 0(M m + M m+1 ) (60) 

in the integrable (= Poisson) case are given by 



7 First examples of integrable discretizations 



(61) 



We use Theorem [I] to investigate the Poisson property of several discretizations. We start 
with the following negative result. 

Proposition 5 The function 7o(M, M, h) = 1 defines a non-Poisson map. 

The simplest way to find Poisson (and therefore integrable) discretizations is to assure that 
the left-hand side of the equation floTf ) does not depend on (p. 

Proposition 6 The functions 

— 4M 2 

7l (M, M, h) = ^ ^- (62) 

(M + M,M + M) V ; 

and 

— - 4E 

72 (M, M, h) = ^ ^— (63) 

(M + M,ri(M + M)) 

define Poisson maps. 
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Proof. Using the second expressions in ( pq )— (|55|) , it is not difficult to derive the following 
formulas: 

™ ™ Ox a 2 + c 2 -6 2 sn 2 (A</?) 

i (M + M,n ( M + M)) = i _^ n2( _ )3n2( ^ ) . 

Hence, the left-hand side of the equation (|50|) for 7 = 7^2 does not depend on <p. Taking 
into account that, according to flTDD, 

a 2 + c 2 = M 2 , ^ + ^ = ^, 
we find the following expressions for this left-hand side: 

1 - ai i2 sn 2 (Ay?) , 

where 

a 1 (M 2 ,E) = ^, a 2 (M 2 ,E) = ^. (64) 

(see ( |10D for the expression of b 2 through M 2 and E). Therefore, the equation floTf ) in these 
two cases is equivalent to 

sn(A^) _hv 1 



cn(Av?)dn(Av9) 2 1 - ai j2 sn 2 (A^) 
Obviously, its solutions satisfy 

A<p = 6 lfl (M*,E,h) = hv{ ^> E) (l + 0(h)). 



8 Veselov— Moser discretization 

In order to introduce the Veselov-Moser discretization, we need, first of all, the matrix 
notation for the Euler top equation (0). Using the well-known isomorphism between the Lie 
algebra (m 3 , x^j with the Lie algebra (so(3), [•, ■]) of 3 x 3 skew-symmetric matrices with the 
usual commutator, we can rewrite the equations of motion (Q) also in the matrix form 

M = [M, Q(M)] , (66) 
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where 



M 



t 

-M 3 

V M 2 



and 



fi(M) 



/ 

-M 3 /C 



M 3 -M 2 
M 1 
-Mt 

M 3 /C -M 2 /£ 
Mi /A 



\ M 2 /B -Mi/ A 
The relation between the matrices M and Q = Q(M) may be expressed as follows: 

M = Jfi + OJ , 

where the entries of the diagonal matrix 

J = diag(Ji, J 2 , J 3 ) 

are defined by the relations 

A = J 2 + J 3 , B = J 3 + Ji, C = J i + J 2 . 
A spectral parameter dependent Lax representation for (|66|) is: 

(M + AJ 2 )' = [M + A J 2 , ft(M) + A J] . 



(67) 



(6J 



(69) 



(70) 



(71) 



(72) 



Now we can describe the Veselov-Moser construction. The differential equation (|66|) is 
replaced by the difference one, 

M = u T Mu , (73) 

where uj G SO (3) is an orthogonal matrix related to M by means of the following relation, 
coming to replace, or to approximate, (|69|): 



hM = uJ - Ju T . 
It is easy to see that the previous two relations imply also 

hM = Ju- lu t J . 



(74) 



(75) 



Moser and Veselov demonstrated that, in general, the equation (f73[) does not determine the 
matrix uj G SO (3) uniquely. This non-uniqueness may be described as follows. Notice that 
(fT4j), (1751) are equivalent to the following factorizations of matrix polynomials: 



I - hXM - A 2 J 2 = (u + XJ){u T - A J) 
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and 

I - ft, AM - A 2 J 2 = (u T - AJ)(w + A J) , 

respectively. Denote by S the set of 6 roots of the equation det(I — hXM — A 2 J 2 ) = 0. 
Then to each decomposition S = S + n 5_ into two disjoint sets of three roots satisfying 
S + = S + = — there corresponds a unique solution uj G SO (3) of (|74| ) such that S+ is the 
set of the roots of det(u; + AJ) = 0, while S_ is the set of the roots of det(u; T — AJ) = 0. 

However, in the continuous limit, when h is supposed to be small, there exists a unique 
decomposition, for which the roots from S_ are positive real numbers 0(h)-close to J-f , 
J 2 _1 , J3" 1 . This is the only decomposition for which the corresponding uj has the following 
asymptotics: 

to = 1 + hQ(M) + 0(h 2 ) . (76) 

Supposing h small enough, we consider from now on only this choice of u), and demonstrate 
that the discretization of Veselov and Moser also belongs to the class of special discretizations 
described by the equation (R9). 

Theorem 2 The Veselov-Moser discretization for h small enough may be represented in the 
form ( |39D with 

7vM (M, M, h) = , 2 == . (77) 

l-i/i 2 ||fi(M + M)|| 2 



Proof. We shall need the following lemma, proof of which is put in the Appendix. 

Lemma 2 For a matrix uj G SO(3), 0(h)-close to I , there exists a unique matrixW G so(3) 
such that W = 0(h), and 



to = I + W + |W 2 

where 

2 

7 " 1 + V 1 - l w l 2 

i/ere 

|W| 2 = Wl + W| + Wl for W = I -W3 Wi I G so(3) . 




With the help of this lemma we derive the following two equations: 

uj - cu T = 2W , uj + uj t = 21 + 7W 2 . (78) 
Now from (ff4|), (ff5|) and the first equation in (|78|) we find: 

/i(M + M) = (w - uj t )J + J(uj - uj t ) = 2(WJ + JW) , (79) 
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hence 



W = -M1(M + M) • 



JO) 



Further, from ([74]), ( [75]) and the second equation in ( f7H| ) we find: 

/i(M-M) = J(cj + a; T ) - (cu + c^J = 7 [J,W 2 ] = 7 [WJ + JW,W] 
= ^ 7 /i 2 [M + M,ft(M + M)] , 



where 



7 



1 



1 - i/i 2 |fi(M + M)p 



x 



and (so (3), 



we see that the 



Using the above-mentioned isomorphism between 
theorem is proved. ■ 

It remains to reproduce from our point of view the result of Moser and Veselov concerning 
the Poisson property of their discretization. To this end, we demonstrate that the equation 
(p0|) with Tvm corresponding to 7vm allows a (^-independent solution for Acp. 

Theorem 3 For the Veselov-Moser discretization the equation (p0|) is equivalent to 



hv 



sn(Ay3)cn(A(/>)dn(Ay3) = — (l + asn 2 (A<p) - /3sn 4 (Ay3) 



(82) 



where 



a 



' 2ACE - (A + C - B)M 2 ' 
(B — C)(AE — M 2 ) 

2ACE- (A + C -B)M 2 



(A — B)(M 2 — CE) 
The solution of ([82]) satisfies 



(3 



AC M 2 -CE ] 



(B - C) 2 AE - M 2 

AC AE - M 2 
k (A - B) 2 M 2 - CE t 



A^ = 5 VM (M 2 , E, h) = Mil f 2;jE;) (l + 0(h)) 



(83) 



(84) 



Proof. For the Veselov-Moser discretization the left-hand side of the equation (|50[) can be 
calculated with the help of (|77j) and the second expressions in (|53|)-(|55l) : 



T(M 2 ,E,p, A<p, h) 



1 Psn 2 (^)sn 2 (A(^) l _ fc 2 sn 2(^) SI1 2(A^) + ^(l _ P sn 2 (^)sn 2 (A^)) 2 - 



(85) 
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where 



a 2 f cn%)cn 2 (A^) 1 b 2 2 2 ^^(A^ + ^ J dn %) dn2 ( A ^) 1 f8m 
G "I^ j dn 2 (^)dn 2 (A^) j + ^ Sn ( ^ )cn (A ^ )dn ^ + ^2 j C n%)cn 2 (Av?) / ' (86) 

Hence the equation (|50|) in the present case is equivalent to: 



1 - Psn%)sn 2 (A^) + Ml - fc 2 sn 2 (^)sn 2 (A^)) 2 - h?G = hv cn ( A ^) dn ( A( ^) 



where 



fc 2 sn 2 (A(/?) ^ 1 ' ' sn 2 (A(/?) 

so that (|87|) is equivalent to 

/i 2 Gi _ cn(Av?)dn(Ay2) 



/c 2 sn 2 (A<£>) sn(Ay9) 
or 

sn(A^)cn(Av9)dn(A^) = . 
This is equivalent to ([32l), as one easily calculates from the definition (|90|) that 

Gx = fcV(l + asn 2 (A<^) -/?sn 4 (A^)" 

18 



sn(Ayj) 

We leave the radical alone on the left-hand side and square the resulting equation to derive: 

- h 2 G = frV ^A^dn 2 ^) _ cn(A y )dn(A y ) , _ ^ X(A } ) (g7) 

sn 2 (Ay9) sn(Ay9) V vr/ v r V v ' 

Obviously, we have: 

G = G {M 2 , E, Ay?) - Gi(M 2 , E, A<p) sn%) , (88) 

a 2 f cn 2 (Ay,) 1 c 2 f dn 2 (Ay,) 1 

Go " ^ 1 dn 2 (A^) / + C~ 2 \ cn 2 (Ay.) j ' (89) 

„ a 2 ( cn 2 (Ao9) 1 c 2 f Fdn 2 (Ao9) 1 b 2 2/ . , , 2/ . . . . 

^ | *»dn*(J£) ) + C? { cn*(A *,) ) " ^(^A^f) ■ 0») 

We shall prove that 

d - A; 2 sn 2 (A^)G = fcVcn 2 (A^)dn 2 (A<^) . (91) 
This formula allows to derive from (jB8|): 

a = - W Wm '(A,)] - ^ s^Wg) , (92) 



with a and (3 as in ([83]). It remains to prove the relation (|9~l"l) . But it follows from the 
definitions 



G x - fcW(Ap)G 



a 2 k 2 c 2 _ b 2 
A 2+ ~C 2 ~~B 2 



k 2 a 2 



> cn 



+ 



I A 2 C 2 B 2 ) 

and the expressions in the curly brackets here in both cases are equal to k 2 i> 2 , due to (|T0|). 



9 Conclusions 

In the present work we have introduced a large family of discretizations of the Euler top 
sharing the integrals of motion with the continuous time system. We characterized those 
of them which are also Poisson with respect to the invariant Poisson bracket of the Euler 
top. For all these Poisson discretizations we found a solution in terms of elliptic functions 
which allows a direct comparison with the continuous time case. We demonstrated that 
the Veselov-Moser discretization also belong to our family, and applied our methods to this 
particular example. 

Let us mention some of the possible dirsections of the further progress. First of all, our 
construction works on the level of reduced equations of motion of the rigid body (in the Lie 
algebra E 3 so (3) ~ su(2)). It would be important to lift it to the level of the cooresponding 
Lie group SO(3) (or SU(2)). On a more mechanical language, we want to discretize the 
equations of motion in the rest frame, and not only in the frame attached firmly to the 
body. Notice that the Veselov-Moser construction has a variational (Lagrangian) origin, 
which makes it valid in the group. Another example of such a discretization admitting 
a variational formulation in the corresponding group is the recently found discrete time 
Lagrange top @. 

As a second problem, our construction has to be generalised to higher dimensions (for the 
general Euler-Manakov case of the multidimensional rigid body). Notice that the Veselov- 
Moser construction goes through in higher dimensions. 

Futher, it would be important to include our construction into the general r-matrix 



framework, connected to factorizations in the loop groups, cf. |S2 . 
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A Proof of Lemma 2 

Let u G SO(3), u = 1 + 0(h). Set V = ^— [, then V G so(3), V = 0(h), and 

1 + V 2 

" = r^v = r^v- 1 - 

But for V G so(3) we have: V 3 = — |V| 2 V, and by induction 

V 2fc+i = (_!)fc|v| 2fe V , V 2k+2 = (-l) fc |V| 2fc V 2 , k > 1 , 

so that 

— = 1 + V+ V 2 , 

I-V 1 + |V| 2 1 + |V| 2 ' 

hence 

^ = I + TTW V + TTW V2 ' 

Now setting 

W = TTW V ' 

so that W G so(3), W = 0(h), we find: 

cj = I + W + -W 2 , where 7 = 1 + |V| 2 . 

We have, obviously, 

2- |W| 2 -2 v /l- |W 



|W| 2 = 0IVI 2 =>■ |V| 2 = 



(1 + |V| 2 ) 21 '' ' ''' |W| 
7 = 1 + |V| 2 



and 

2-2^1 - |W 



W 



l + ^/l- |W| 



This proves the lemma. 
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